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Abstract 

Fuel  cells  are  being  considered  as  an  important  technology  that  can  be  used  for  various  power  applications.  For  portable  electronic  devices 
such  as  laptops,  digital  cameras,  cell  phone,  etc.,  the  direct  methanol  fuel  cell  (DMFC)  is  a  very  promising  candidate  as  a  power  source. 
Compared  with  conventional  batteries,  DMFC  can  provide  a  higher  power  density  with  a  long-lasting  life  and  recharging  which  is  almost 
instant.  However,  many  issues  related  to  the  design,  fabrication  and  operation  of  miniaturized  DMFC  power  systems  still  remain  unsolved. 
Fuel  delivery  is  one  of  the  key  issues  that  will  determine  the  performance  of  the  DMFC.  To  maintain  a  desired  performance,  an  efficient 
fuel  delivery  system  is  required  to  provide  an  adequate  amount  of  fuel  for  consumption  and  remove  carbon  dioxide  generated  from  fuel  cell 
devices  at  the  same  time.  In  this  paper,  a  novel  fuel  delivery  system  combined  with  a  miniaturized  DMFC  is  presented.  The  core  component 
of  this  system  is  a  piezoelectric  valveless  micropump  that  can  convert  the  reciprocating  movement  of  a  diaphragm  activated  by  a  piezoelectric 
actuator  into  a  pumping  effect.  Nozzle/diffuser  elements  are  used  to  direct  the  flow  from  inlet  to  outlet.  As  for  DMFC  devices,  the  micropump 
system  needs  to  meet  some  specific  requirements:  low  energy  consumption  but  a  sufficient  fuel  flow  rate.  Based  on  theoretical  analysis,  the 
effect  of  piezoelectric  materials  properties,  driving  voltage,  driving  frequency,  nozzle/diffuser  dimension,  and  other  factors  on  the  performance 
of  the  whole  fuel  cell  system  will  be  discussed.  As  a  result,  a  viable  design  of  a  micropump  system  for  fuel  delivery  can  be  achieved  and  some 
simulation  results  will  be  presented  as  well. 
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1.  Introduction 

Many  portable  electronic  devices  available  to  consumers 
nowadays  are  combining  more  and  more  functionalities  to¬ 
gether.  Consequently,  the  power  consumption  of  these  de¬ 
vices  is  consistently  increasing.  Although  a  great  amount 
of  effort  has  already  been  put  into  the  research  on  the  im¬ 
provements  of  conventional  rechargeable  batteries,  the  per¬ 
formance  of  these  power  devices  is  still  unsatisfactory.  On 
the  other  hand,  as  reported  by  Hebling  et  al.  [1],  portable  and 
miniaturized  fuel  cells  have  shown  many  advantages  over  the 
batteries,  such  as  higher  power  density,  instant  recharging, 
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low  environmental  impacts  etc.  So,  it  becomes  very  possible 
to  use  fuel  cell  as  a  portable  power  source  instead  of  batteries 
in  the  near  future. 

Among  various  types  of  fuel  cells,  the  proton  exchange 
membrane  fuel  cell  (PEMFC)  and  DMFC  are  applicable  for 
the  portable  electronics  due  to  their  low  operation  tempera¬ 
ture.  In  comparison  with  DMFC,  the  PEMFC  device  that  uses 
hydrogen  directly  as  fuel  has  a  better  performance,  but  suf¬ 
fers  from  the  storage  problem  of  hydrogen.  So  far  the  storage 
and  transportation  of  hydrogen  are  still  not  cost-effective  for 
commercial  applications.  And  the  fuel  cell  with  a  reformer  to 
convert  the  liquid  fuel  to  hydrogen  gas  is  also  less  favorable 
because  it  increases  the  complexity  of  the  system  integration 
and  generates  undesired  byproduct  (carbon  monoxide,  CO) 
during  the  reforming  process.  Therefore,  a  DMFC  using  a  liq- 
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uid  methanol  solution  as  fuel  directly  is  much  more  straight¬ 
forward  and  attractive.  Many  researchers  including  Waidhas 
et  al.  [2]  and  Scott  et  al.  [3]  have  discussed  the  benefits  and 
challenges,  utilizing  such  a  direct-fuelled  system. 

To  optimize  a  miniaturized  DMFC  power  system,  many 
engineering  issues  need  to  be  considered  which  include  ther¬ 
mal  management,  methanol  crossover,  fuel  delivery,  oxidant 
supply  and  power  conditioning.  This  paper  mainly  focuses  on 
the  fuel  delivery.  For  medium-  and  large-size  DMFC,  con¬ 
ventional  pumps  can  be  used  to  drive  the  fuel  flowing  through 
the  fuel  cell.  For  miniaturized  DMFC,  the  flow  rate  of  fuel 
is  very  small  and  there  are  several  options  for  fuel  delivery. 
Mench  et  al.  [4]  presented  a  pump-less  and  self-activated  de¬ 
sign  for  a  DMFC  stack,  in  which  gravitational  and  capillary 
forces  were  used  to  feed  the  liquid  methanol  solution.  By 
their  calculation,  the  fuel  supply  rate  was  adequate  for  ac¬ 
ceptable  performance.  This  design  was  quite  interesting  and 
simple  but  no  experimental  results  were  reported.  In  addition, 
the  fabrication  of  the  DMFC  stack  seems  to  be  very  compli¬ 
cated.  The  more  general  way  is  to  choose  a  micropump.  Both 
Shukla  et  al.  [5]  and  Scott  et  al.  [6]  used  a  peristaltic  pump 
in  their  experimental  study.  The  flow  rate  of  fuel  varied  from 
less  than  1  to  5  ml  min-1  and  the  active  cell  area  ranged  from 
9  to  25  cm2.  The  operation  performance  of  their  fuel  cell  sys¬ 
tems  was  in  the  acceptable  range.  Therefore,  it  seems  to  be 
more  practical  and  easier  to  employ  a  micropump  for  fuel 
delivery. 

The  research  on  micropumps  started  from  1980s  and  var¬ 
ious  designs  have  been  presented  since  then.  Olsson  [7] 
summarized  the  performance  of  several  typical  micropumps. 
Among  them  the  reciprocating  diaphragm  pumps  show  a 
broad  range  of  pumping  flow  rate  and  can  be  easily  fabri¬ 
cated  using  a  MEMS  technique.  To  eliminate  the  wearing 
and  fatigue  of  valves  and  reduce  clogging,  a  valveless  de¬ 
sign  is  preferred.  Stemme  E.  and  Stemme  G.  [8]  introduced 
the  first  piezoelectric  valveless  micropump.  Gerlach  et  al.  [9] 
presented  a  similar  type  of  micropump  using  different  flow 
directing  elements.  Forster  et  al.  [10]  compared  the  perfor¬ 
mance  between  the  micropump  with  diffuser  and  the  microp¬ 
ump  with  valvular  conduit.  It  was  found  that  the  valvular  con¬ 
duit  design  had  a  higher  volumetric  efficiency.  Olsson  [7]  and 


Ullmann  [11]  also  conducted  fairly  extensive  modeling  and 
experimental  research  on  this  valveless  micropump.  Based  on 
these  previous  studies,  a  new  fuel  delivery  system  for  minia¬ 
turized  DMFC  using  a  piezoelectric  valveless  micropump  is 
presented  in  this  paper.  The  performance  of  this  system  is 
analyzed  theoretically. 

2.  System  design 

The  cross-section  of  the  DMFC  power  system  with  piezo¬ 
electric  micropump  is  shown  schematically  in  Fig.  1.  Note 
that  the  figure  is  not  drawn  in  the  actual  scale.  This  system 
mainly  consists  of  the  following  parts:  fuel  cell  membrane 
electrode  assembly  (ME A),  fuel  chamber,  nozzle/diffuser, 
micropump  and  pump  chamber,  and  fuel  supply  manifold. 
All  these  parts  are  fabricated  in  a  multilayer  structure  to  ob¬ 
tain  a  compact  fuel  cell  system.  The  fuel  cell  MEA  is  made 
of  a  Nafion®  117  membrane  layer  sandwiched  by  two  elec¬ 
trode  layers  with  catalysts  deposited  on  them.  And  the  mi¬ 
cropump  is  fabricated  by  adhering  a  thin  piezoelectric  disk 
on  to  a  metal  diaphragm.  When  applying  an  alternating  volt¬ 
age  to  the  piezoelectric  disk,  the  diaphragm  is  actuated  to 
produce  bending  deformation  that  causes  volume  change  in 
pump  chamber.  By  selecting  an  appropriate  shape  and  di¬ 
mension  of  the  nozzle/diffuser  between  the  pump  chamber 
and  the  fuel  chamber,  the  fuel  can  circulate  in  the  desired 
direction.  And  the  fuel  supply  from  the  right  chamber  can 
compensate  for  the  fuel  consumption. 

3.  Model  of  micropump 

As  shown  in  Fig.  2,  the  vibrating  diaphragm  on  top  of 
the  pump  chamber  is  actually  a  piezoelectric  bending  mode 
actuator  with  fixed  edges.  Both  applied  voltage  across  the 
PZT  disk,  U,  and  pressure  difference  exerted  on  the  whole 
diaphragm,  p ,  can  cause  the  deflection  in  vertical  direction. 
Specifically  p  is  the  difference  between  the  pressure  above 
the  pump  diaphragm  (usually  the  atmospheric  pressure)  and 
the  pressure  below  the  pump  diaphragm  (the  pressure  of  the 
pump  chamber). 
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Fig.  1.  Schematic  of  the  miniaturized  DMFC  system  driven  by  piezoelectric  valveless  micropump. 
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Fig.  2.  Schematic  of  the  piezoelectric  micropump  actuator. 


A  theoretical  model  developed  by  Li  and  Chen  [12]  is 
adopted  to  estimate  the  deflection  induced  by  the  applied 
voltage.  According  to  their  study,  this  analytical  solution  is 
able  to  predict  reasonable  deflections  for  the  passive  plates 
of  various  Young’s  modulus  from  30  to  600  GPa.  The  passive 
plate  used  in  this  design  is  well  within  this  range.  To  simplify 
the  modeling  process,  the  piezoelectric  pump  actuator  is  di¬ 
vided  into  two  parts:  one  is  the  three-layer  structure  including 
the  PZT  disk,  the  bonding  layer  and  the  passive  metal  plate; 
the  other  is  the  rest  part  of  the  passive  plate. 

From  Li  and  Chen  [12],  the  deflections  of  the  three-layer 
structure  and  the  rest  part  of  the  passive  plate  can  be  calculated 
separately.  For  the  rest  part  of  the  passive  plate,  its  deflection 
W\  is  obtained  using  the  equations  given  by  Timoshenko  et 
al.  [13] 


W\  (r)  = 


2Dp[(l 


Mnb2 

Yp)ci2  +  (1  +  Yp)b2] 


X 


+  2 a2  log  — |-  a 
a 


For  the  three-layer  structure,  its  deflection  W2  also  follows 
Timoshenko  et  al.  [13] 


Mnb2 

W?(r)  = - — - 

2DP[(1  -  YP)a2  +  (1  +  YP)b2] 

x  (—b2+2a2  log  -+a2  \  +  M'  -ib2-r2) 

\  a  )  2Dc(1+Yc) 

(2) 

Here,  a  and  b  represent  the  radius  of  the  passive  plate  and  the 
PZT  disk  respectively.  Dp  and  yp  are  the  flexural  modulus  and 

Poisson’s  ratio  of  the  passive  plate.  Dc  and  yc  are  the  flexural 
modulus  and  Poisson’s  ratio  of  the  three-layer  structure.  M\ 
and  Mu  are  intermediate  bending  moments. 


The  bending  moment  Mi ,  as  derived  by  Li  and  Chen  [12], 
is 


r]Ud^i(2h'  +  2h\y  +  hp7x)(hr^  +  h 3) 

2 [h3  +  (1  -  a)h ,3  +  (a  -  p)(ti  +  hb )3 

+/3(h'  +  hb  +  /*pzt)3] 


Here,  U  is  the  static  voltage  applied  on  the  PZT  disk,  d^\  is 
the  piezoelectric  charge  coefficient,  and  h  is  the  location  of 
the  neural  surface  of  the  three-layer  structure.  The  equations 
for  other  variables  are: 


a  = 


1  -  Xp  Eb 

1-Kb2  V 


1  Vp  £Pzt 

1  -  Vp/.i  Ep 


Epzt 
1  —  Vpzt 


where  Ep,  E^,  and  Epzt  represent  the  Young’s  modulus  of  the 
passive  plate,  the  bonding  layer  and  the  PZT  disk,  respec¬ 
tively.  y p,  y b,  and  ypzt  represent  the  Poisson’s  ratio  of  the 
passive  plate,  the  bonding  layer  and  the  PZT  disk,  respec¬ 
tively.  hp,  hb,  and  hpzt  represent  the  thickness  of  the  passive 
plate,  the  bonding  layer  and  the  PZT  disk,  respectively. 

Based  on  the  work  done  by  Dobrucki  et  al.  [14],  the  loca¬ 
tion  of  the  neutral  surface  of  the  three-layer  structure,  h,  can 
be  calculated  by 


Eph2/(\  -  Kp2)  +  Eb[(hp  +  hb)2  -  h2]/(l  -  yI) 

1  "I-  £pzt[(6P  /fpzt)  (bp  +  hb)  ]/( 1  Ypzt) 

2  Ephp/(l-Y$)+Ebhb/(l-Yb)+EpZthpZt/(l-YpZt) 

(4) 


By  observing  Eqs.  (1)  and  (2),  it  can  be  found  that  the 
continuity  of  deflection  at  r  =  b  is  already  satisfied.  So,  the 
continuity  condition  of  the  slope  of  the  deflection  at  r-b  has 
to  be  used  to  determine  unknown  moment  M\  \ . 
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To  satisfy  the  continuity  condition, 


dW\(r) 

dr 


r—b 


dW2  (r) 
dr 


r=b 


Substituting  Eqs.  (1)  and  (2)  into  Eq.  (5),  the  moment  M\\  is 
found  to  be 


Mn 


MiDp[(l  -  yp)a2  +  (1  -  yp)b 2] 

Dc(l  +  Yc)(a2  ~  b 2) 


Substituting  Eq.  (6)  back  into  Eqs.  (1)  and  (2),  the  whole 
deflection  can  be  rewritten  as 


Wi(r)  = 


M\b2{r2  —  2  a2  1  ogr/a  —  a2) 
2DC(1  +  yc)(a2  ~  b2) 


(b  <  r  <  a)  (7) 


M\  [(az —bz)(bz —rz)+b2(b2—2az  log  b/a-az)] 


W2(r)= 


2Dc(\+yc)(a2  —  b2) 


(0  <  r  <  b) 


(8) 


Fig.  3.  Comparison  of  the  diaphragm  deflections  under  100  V  and  1  kPa. 


The  deflection  caused  by  pressure  difference  can  be  eval¬ 
uated  using  the  similar  method.  Under  constant  pressure  dif¬ 
ference  p ,  the  deflection  of  the  rest  part  of  the  passive  plate 
is 


W3(r)  = 


P  (a2  —  r2)2 


64  D 


P 


(M21  —  M2)b2(r2  —  2  az 1  ogr/a  —  az) 


+ 


2£>p[(l  -  yv)a2  +  (1  +  yv)b2] 
The  deflection  of  the  three-layer  structure  is 


W4(r)  = 


P  ,2  7  2'  2 


64  D 


(a  —  b  ) 


P 


(M21  —  M2)b2(b2  —  2  az  1  ogb/a  —  az) 
2Dp[(l  -  yv)a2  +  (1  +  yv)b 2] 

P{h2  ~  F)  /  5  +  Ycb2  _  r: 

64Z)p  \  1  +  yc 

Ml  ( b 2  -  r2) 


+ 


+ 


+ 


2Dc(l  +  yc) 


(9) 


(10) 


Here,  M2  and  M21  are  also  two  intermediate  bending  mo¬ 
ments,  of  which 


M21  =  ^[(!  +  Yp)a2  -  (3  +  Yp)b 2] 


(11) 


Similar  to  the  case  of  the  deflection  due  to  piezoelectric 
effect,  the  continuity  of  deflection  at  r  =  b  is  already  satisfied 
in  Eqs.  (9)  and  (10).  And  the  continuity  condition  of  the  slope 
of  the  deflection  at  r-b  can  be  used  to  determine  the  other 
moment  M2. 

The  volume  change  of  the  pump  chamber  induced  by  the 
deflection  of  the  pump  diaphragm  can  be  calculated  using 


AE  =  2  ix 


R 

W(r)r  dr 

0 


(12) 


The  total  deflection  of  the  pump  diaphragm  is  the  superpo¬ 
sition  of  the  deflection  introduced  by  applied  voltage  and  the 
deflection  caused  by  pressure  difference.  In  actual  pump  ap¬ 
plications,  the  applied  voltage  can  be  easily  controlled  while 
the  pressure  of  the  pump  chamber  is  coupled  with  the  deflec¬ 
tion  of  the  pump  diaphragm  and  flow  resistance.  Numerical 
iteration  methods  can  be  used  to  determine  the  deflection  of 
the  pump  diaphragm  precisely. 

For  this  design,  several  kPa  of  pressure  difference  between 
the  fuel  chamber  and  the  pump  chamber  should  be  enough 
to  drive  the  fuel  circulating  between  them.  If  the  pressure 
of  the  fuel  chamber  is  kept  at  the  atmospheric  pressure,  the 
pressure  difference  between  the  pump  chamber  and  the  at¬ 
mospheric  pressure  is  also  in  the  range  of  several  kPa.  Using 
Eqs.  (7)-(10),  the  deflection  under  100  V  of  applied  voltage 
and  the  deflection  under  1  kPa  of  pressure  difference  are  cal¬ 
culated  separately  and  the  results  are  shown  in  Fig.  3.  The  di¬ 
mensions  and  material  properties  of  the  pump  actuator  listed 
in  Table  1  are  used  for  calculation.  In  comparison  with  the  de¬ 
flection  induced  by  applied  voltage,  the  deflection  caused  by 
pressure  difference  is  small  and  the  resultant  volume  change 
is  smaller.  Therefore,  the  effect  of  the  pressure  difference  on 
the  volume  change  is  neglected  to  obtain  an  analytical  solu¬ 
tion. 


Table  1 


Dimensions  and  material  properties  of  the  piezoelectric  micropump  actuator 


Piezoelectric 
disk  (PZT-5H) 

Bonding  layer 
(conductive  epoxy) 

Passive 
plate  (steel) 

Diameter  (mm) 

20 

20 

24 

Thickness  (mm) 

0.2 

0.02 

0.25 

Young’s  modulus 

6.2  X  1010 

5.17  x  109 

1.95  X  1011 

(Pa) 

Poisson’s  ratio 

0.31 

0.3 

0.3 

d3i  (mV-1) 

-3.2  X  10“10 

£^3  (Fm-1) 

2.856  x  10“8 
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Due  to  this  simplification,  the  total  deflection  of  the  pump 
diaphragm  can  be  expressed  as 


the  same  flow  velocity.  The  characteristic  of  nozzle/diffuser 
element  can  be  described  as  follows, 


W(r)  = 


W2(r),  0  <  r  <  b; 

and  R  =  a. 

W\(r),  b  <  r  <  a; 


(13) 


Substituting  Eq.  (13)  into  Eq.  (12),  the  volume  change  of  the 
pump  chamber  is  obtained 


AV  = 


TC 

2 


kpJjA  +  2k^b2  +  kc 


b 4  +  4 a2b2  log 


(14) 


Here, 


M\  M\b2(b2  —  2 a2  log  b/a  —  a2) 

2DC(1  +  yc)  ’  B  =  2DC(1  +  yc)(a2  -  b2) 

M\b2 

2DC{\  +  yc)(a2  -  b2) 


It  is  obvious  that  these  coefficients  k\,  ,  kc ,  are  proportional 
to  M i .  On  the  other  hand,  M\  is  proportional  to  the  applied 
voltage  U.  Therefore,  the  volume  change  AV  is  ultimately 
proportional  to  the  applied  voltage  U,  which  can  be  written 
as 


AV  =  kQffU 


(15) 


The  proportional  coefficient  is  given  by 


^eff  — 


TC 

2kv 


k'Ab4+2k'Bb2+k'c  (  a4 


b4  +  4 a2b2  log  - 

a 


(16) 


Here, 


1  ,  b2(b2  —  2 a2  log  b/a  —  a2) 

2Z)C(1  +  yc)  ’  B  =  2DC(1  +  yc)(«2  -  W) 

b 2 

2DC(\  +  Yc)(a2  -  b2)’ 


,  r,d3i(2h' +  2hb  +  hpzt)(h'3  +  h3) 

k\j  = - - - i - - - r 

2[h3  +  (1  -  a)h'3  +  (a  -  P)(ti  +  hby 

+  P(h'  +  hb  +  /zpZt)3] 

If  an  alternating  sinusoidal  voltage  is  applied  to  the  piezo¬ 
electric  disk  and  the  actuation  frequency  is  far  lower  than 
the  resonance  frequency  of  the  pump,  U  can  be  substituted 
by  Um  sin  (cot).  Um  is  the  amplitude  of  the  sinusoidal  volt¬ 
age  applied  on  the  piezoelectric  plate  and  co  is  the  angular 
frequency. 

To  estimate  the  net  flow  rate  of  the  pump,  the  modeling  of 
the  nozzle/diffuser  is  also  necessary.  As  mentioned  before, 
these  nozzle/diffuser  elements  can  directly  flow  from  inlet 
to  outlet.  They  are  geometrically  designed  to  have  a  lower 
pressure  loss  in  one  direction  than  in  the  opposite  direction  for 


A  p  =  \pv2%  (17) 

Here,  A p  is  the  pressure  loss  through  the  nozzle/diffuser  el¬ 
ement,  p  is  the  density  of  the  fluid,  v  is  the  velocity  of  the 
fluid,  and  §  is  the  pressure  loss  coefficient.  And  it  can  also  be 
expressed  in  the  form  of  flow  rate  Q , 


Q  =  CfAp 


(18) 


Here,  C  is  called  conductivity  coefficient.  Two  different  flow 
directions  correspond  to  two  different  C  values.  One  is  higher 
than  the  other.  The  conductivity  coefficient  in  positive  direc¬ 
tion  is  represented  by  Ch  and  the  conductivity  coefficient  in 
negative  direction  is  represented  by  Cl.  Accordingly,  there 
are  also  two  different  §  values:  the  pressure  loss  coefficient 
in  positive  direction  is  represented  by  ^positive  and  the  con¬ 
ductivity  coefficient  in  negative  direction  is  represented  by 

^negative- 

As  mentioned  before,  it  is  assumed  that  the  effect  of  the 
pressure  difference  on  the  volume  change  is  negligible,  and 
thus  only  the  volume  change  induced  by  applied  voltage  is 
considered.  Also  assume  that  there  is  no  pressure  variance 
throughout  each  chamber  so  that  the  pressures  of  the  pump 
chamber  and  the  fuel  chamber  can  be  represented  as  pv  and pf , 
respectively.  Given  these  assumptions,  the  equations  derived 
by  Ullmann  [15]  can  be  used  to  obtain  the  analytical  solutions 
for  the  fuel  flow  rate  and  the  pressure  of  the  pump  chamber. 
The  pressure  of  the  fuel  chamber  is  set  to  be  the  atmospheric 
pressure. 

(1)  Pump  mode:  pv  >  pf 


01  = 


02  = 


CL 


Ch  +  Cl 


Ch 


Ch  +  Cl 


keff Umco  cos  (cot) 


kQffUmco  cos  (cot) 


(19) 


(20) 


Pv  =  Pf  + 


kQffUm<z>  cos  (cot) 
Ch  +  Cl 


i  2 


(21) 


(2)  Supply  mode:  pv  <  pf 


01  =  — 


02  =  — 


Ch 


Ch  +  Cl 


Cl 


Ch  +  Cl 


keffUm00  cos  (cot) 


kQffUmco  cohort) 


(22) 


(23) 


Pv  =  Pf 


keff  U m (O  COS((Ut) 

Ch  +  Cl 


i  2 


(24) 


The  average  net  flow  rate  is  given  by 


0  = 


ke,ffUmco  Ch  Cl  kQfvUmco  t]  ^2  1 


TC  Ch  +  Cl 


TC  Y)1!2  +  1 


(25) 


r 

Here,  CJ  =  ^negative/ ^positive  =  (Ch/ Cl) 
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Since  the  electrical  power  used  to  drive  the  micropump 
comes  from  the  fuel  cell  itself,  it  is  necessary  to  estimate  how 
much  energy  the  micropump  needs  to  consume.  Electrically 
the  bending  piezoelectric  actuator  behaves  like  a  planar  ca¬ 
pacitor.  The  electric  capacitance  under  the  given  mechanical 
boundary  condition  is  very  complicated.  In  a  similar  case,  the 
capacitance  of  a  circular  PZT  actuator  with  clamped  edges 
has  given  by  Kim  et  al.  [16] 

cT  jrh 2 

q  _  £33nD 
hpzt 


(26) 


3(^n)  Sp^p^pztC^pzt+^p) 
ShB  31 


1- 


K 


31 


Here,  is  the  permittivity  of  the  PZT  disk.  is  the  elastic 
compliance  of  the  PZT  disk  at  constant  electric  field  and 
sv  is  the  elastic  compliance  of  the  passive  layer.  K^i  is  the 
electromechanical  coupling  coefficient  of  the  PZT  disk,  v  is 
the  Poisson’s  ratio.  Kim  et  al.  [16]  assume  that  the  Poisson’s 
ratio  of  the  PZT  disk  is  the  same  as  that  of  the  passive  plate  and 
the  effect  of  the  bonding  layer  on  the  capacitance  is  neglected. 
The  equations  for  the  other  two  parameters  are: 

6*h  =  ^pzUp  T  ^pUl’ 

^31  ^pzUp  "3"  4snsvhvhpzt  +  6^11^p^p^pzt 
+4^1Sp/jp/?pzt  +  /?pOn)2 


The  term  outside  the  parenthesis  of  Eq.  (26)  is  actually  the 
capacitance  with  free  boundary  condition.  Using  the  dimen¬ 
sions  and  material  properties  of  the  piezoelectric  actuator  in 
our  design  for  calculation,  it  is  found  that  the  capacitance  with 
clamped  boundary  condition  is  about  30%  less  than  the  capac¬ 
itance  with  free  boundary  condition.  Although  the  boundary 
condition  of  our  case  is  not  exactly  the  same  as  that  of  Kim  et 
al.  [16],  the  results  of  two  cases  should  be  close.  For  approx¬ 
imation  the  capacitance  with  free  boundary  condition  is  used 
here  because  the  resultant  estimation  of  the  pump  power  con¬ 
sumption  will  be  larger,  which  is  safe  for  our  system  design. 
The  capacitance  with  free  boundary  condition  is  described  as 


permittivity  x  area 
distance 


hpzt 


(27) 


Suppose  the  applied  voltage  is  sinusoidal,  that  is 


U  ( t )  =  Um  sin  (cot) 

The  current  is  then, 
d  U(t) 


m  =  C0 


d  t 


(28) 


(29) 


So,  the  instantaneous  power  consumption  is  given  by 
Welec  =  i(t)U(t)  =  ±®C0t/2  sin(2fitf)  (30) 


It  should  be  noted  that  as  a  capacitor,  the  piezoelectric  actua¬ 
tor  would  store  a  large  portion  of  the  input  electrical  energy; 
therefore,  only  part  of  the  input  energy  can  be  converted  into 
output  mechanical  energy  for  fuel  delivery.  The  stored  en¬ 
ergy  will  remain  in  the  driving  circuit  and  will  be  used  in 
next  driving  cycle.  For  the  worst  case,  the  maximum  value 
of  Weiec  is  chosen  to  estimate  the  efficiency  of  the  whole  fuel 
cell  system.  Combining  Eq.  (25)  and  Eq.  (30)  by  eliminating 
Um,  we  have 


max  Weiec  = 


TrCp  f  r]l/2  +  1 

4*4/  V1/2  -  1 


Here, /is  the  frequency  of  the  driving  voltage. 


(31) 


4.  Performance  analysis  of  DMFC 

To  analyze  the  performance  of  DMFC  with  a  micropump 
fuel  delivery  system,  a  model  describing  the  relationship  be¬ 
tween  the  flow  rate  of  the  liquid  fuel  and  fuel  cell  J-V  curve 
would  be  more  convenient.  However,  for  simplicity,  we  will 
use  some  appropriate  experimental  data  from  literature  to 
discuss  the  performance  of  fuel  cells.  Up  to  now  some  ex¬ 
perimental  data  have  been  reported  about  the  performance 
of  DMFC  under  different  operation  conditions.  For  example, 
Jung  et  al.  [17]  evaluated  a  single  cell  (7.6  cm2)  performance 
over  the  range  of  60-120  °C  with  methanol  concentrations  of 
0. 5-4.0  M,  flow  rates  of  9  ml  min- 1 ,  and  an  oxygen  saturation 
of  105  seem.  A  current  density  of  230  mA  cm-2  at  0.55  V  was 
obtained  from  2.5  M  methanol.  Baldauf  et  al.  [18]  conducted 
their  experiments  under  air  operation  at  over  low  pressure. 
A  current  density  of  100  mA  cm-2  and  a  0.5  V  voltage  at 
80  °C  were  observed  from  0.5  M  methanol.  The  electrode 
area  of  their  experimental  cells  was  3  cm2.  Dohle  et  al.  [19] 
reported  an  optimum  methanol  concentration  of  2M  for  a 
20  cm2  single  cell  under  1 10  °C.  The  current  density  and  the 
corresponding  cell  voltage  were  in  a  range  similar  to  that  was 
used  by  Baldauf  et  al.  [18].  The  experiments  performed  by 
Amphlett  et  al.  [20]  indicated  that  there  is  a  complex  relation 
between  the  effect  of  flow  velocity  and  the  influence  of  the 
production  rate  of  CO2. 

The  latest  experimental  data  on  DMFC  were  presented  by 
Gurau  et  al.  [21].  They  measured  a  series  of  J-V  character¬ 
istic  curves  for  a  single  fuel  cell  assembly  (5  cm2)  operating 
under  different  operation  temperatures  (80,  60  and  40  °C), 
different  methanol  flow  rates  (5,  0.5  and  0.15  ml  min-1),  and 
different  methanol  concentrations  (2,  1  and  0.5  M).  The  per¬ 
formance  analysis  in  this  paper  is  mainly  based  on  this  group 
of  data.  From  the  data  they  obtained,  it  is  found  that  at  a 
low  current  density  (less  than  100  mA  cm-2),  the  influence 
of  changing  methanol  concentration  (from  0.5  to  2  M)  on  the 
characteristic  J-V  curve  of  DMFC  is  not  significant.  The  same 
phenomenon  is  also  observed  from  the  experimental  data  ob¬ 
tained  by  other  researchers.  Therefore,  in  our  new  design, 
initially  a  2  M  methanol  solution  can  be  added  to  the  cham- 
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Fig.  4.  Selected  experimental  J-V  data  of  a  single  fuel  cell  at  different  fuel 
flow  rate  (Gurau  et  al.  [21]). 

bers.  During  operation,  the  concentration  of  the  methanol 
solution  will  decrease  due  to  the  fuel  consumption.  But  the 
performance  of  the  fuel  cell  will  not  drop  too  much  if  the  cur¬ 
rent  density  is  not  very  high  and  the  methanol  concentration 
is  still  larger  than  0.5  M.  A  certain  range  of  methanol  concen¬ 
tration  can  be  maintained  through  continuous  fuel  flow  from 
the  fuel  supply  manifold. 

5.  Results  and  discussions 

Fig.  4  shows  the  experimental  data  of  a  single  fuel  cell 
running  at  40  °C.  The  concentration  of  methanol  solution  is 
1  M.  The  specific  data  are  chosen  for  calculation  because  the 
operation  temperature  is  very  close  to  that  of  the  commercial 
applications  and  the  concentration  is  in  the  middle  region. 
Based  on  Fig.  4,  the  power  density  can  be  calculated  as  shown 
in  Fig.  5.  The  maximum  value  is  reached  when  the  current 
density  is  around  120  mA  cm-2  and  the  cell  voltage  is  around 


Fig.  5.  Calculated  power  density  of  a  single  fuel  cell  at  different  fuel  flow 
rate. 


0.35  V.  If  our  design  is  set  to  be  running  under  this  condition, 
for  a  1  W  fuel  cell  assembly,  the  active  area  needed  is  then 
25  cm2  (5  cm  x  5  cm)  and  the  total  current  is  about  3  A. 

The  next  thing  we  need  to  know  is  the  least  fuel  flow  rate 
required  under  this  operation  condition.  As  is  well  known, 
the  electrochemical  reactions  of  a  DMFC  are: 

anode  reaction:  CH3OH  +  H2O  — >►  CO2  +  6H+  +  6e_ 

(32) 

cathode  reaction:  5O2  +  6H+  +  6e-  — >►  3H2O  (33) 

overall  reaction:  CH3OH  +  |02  — >  CO2  +  2H2O  (34) 

As  shown  in  the  above  equations,  ideally  1  mol  methanol  will 
react  with  1  mol  of  water  and  generate  6  mols  of  electrons. 
Therefore,  the  relation  between  the  total  generated  current 
and  the  flow  rate  of  consumed  methanol  solution  can  be  de¬ 
scribed  as 

I  =  Q  x  Dm  x  6  x  F  (35) 

Here,  Dm  is  the  mol  density  of  the  methanol  solution  and  F 
is  Faraday’s  constant  (96485  C  mol-1).  Based  on  Eq.  (35), 
for  1  M  methanol,  the  least  fuel  flow  rate  required  is  approx¬ 
imately  0.3 1  ml  min-1 . 

The  material  properties  and  the  dimensions  of  the  mi¬ 
cropump  actuator  are  listed  in  Table  1.  As  mentioned  be¬ 
fore,  there  are  two  different  designs  of  nozzle/diffuser  ele¬ 
ments  that  can  be  used  to  direct  the  flow.  Stemme  E.  and 
Stemme  G.  [8]  have  used  one  in  their  valveless  micropump, 
of  which  the  opening  angles  are  small,  normally  less  than  20° 
and  the  diffuser  direction  is  the  positive  flow  direction.  Ger- 
lach  et  al.  [9]  fabricated  the  other  in  silicon  using  anisotropic 
wet  etching.  The  opening  angle  of  this  design  is  70.5°  and 
the  positive  flow  direction  is  in  the  converging  wall  direc¬ 
tion.  In  our  design,  a  conical  nozzle  element  with  sharp  in¬ 
let  and  outlet  is  chosen  due  to  its  easier  fabrication  process, 
although  the  performance  is  not  the  best.  The  opening  an¬ 
gle  of  this  nozzle  is  70.5°.  The  minimum  diameter  is  1  mm 
and  the  maximum  diameter  is  3.15  mm.  By  following  the 
same  analysis  method  presented  by  Olsson  [7],  ^positive  and 
^negative  of  this  nozzle  are  found  to  be  1 .005  and  1 .46,  respec¬ 
tively.  The  ratio  of  the  pressure  loss  coefficients,  77,  is  then 
about  1.45. 

Since  the  dimensions  and  materials  of  the  pump  actua¬ 
tor  and  nozzle/diffuser  elements  have  been  determined,  the 
actuation  voltage  needed  to  meet  the  specific  flow  rate  re¬ 
quirement  can  be  calculated  by  Eq.  (25).  For  example,  given 
the  actuation  frequency  is  100  Hz  and  the  fuel  flow  rate  is 
lml  min-1,  Um  is  then  equal  to  71V.  By  using  Eq.  (31), 
the  maximum  electrical  power  consumption  is  about  70  mW, 
which  is  1%  of  the  total  power  produced  by  the  fuel  cell. 
Comparing  with  a  similar  micropump  developed  by  Bardell 
et  al.  [22]  that  can  provide  a  water  flow  rate  up  to  1  ml  min-1 , 
while  consuming  approximately  50  mW  power,  this  estima¬ 
tion  is  not  far  away  from  the  real  case.  Furthermore,  some 


T.  Zhang,  Q.-M.  Wang  /  Journal  of  Power  Sources  140  (2005)  72-80 


79 


Current  density  (mA/cm2) 

Fig.  6.  Net  power  output  under  different  operation  conditions  (actuation 
frequency  =  100  Hz). 

methods  that  will  be  discussed  later  can  be  applied  to  de¬ 
crease  the  power  consumption  of  the  micropump. 

Fig.  6  shows  the  net  output  power  by  subtracting  the  power 
consumed  by  the  piezoelectric  pump  from  the  total  power 
generated  by  fuel  cell  under  different  operation  conditions  at 
a  driving  frequency  of  100  Hz.  Data  from  Fig.  4,  which  was 
obtained  from  the  experiments  on  a  5  cm2  fuel  cell,  are  used 
for  calculations.  So,  the  results  may  not  be  precise  but  the 
trend  should  be  similar.  The  negative  output  power  means  that 
the  power  generated  by  fuel  cell  is  not  enough  to  drive  the  mi¬ 
cropump.  It  is  obvious  that  the  fuel  flow  rate  of  0.5  ml  min-1 
is  preferred. 

From  Eq.  (31),  it  can  be  found  that  increasing  the  driv¬ 
ing  frequency  while  maintaining  the  same  flow  rate  can  re¬ 
duce  the  electrical  power  consumption  of  the  piezoelectric 
micropump.  Fig.  7  shows  the  output  power  curves  at  a  driv¬ 
ing  frequency  of  200  Hz.  The  fuel  cell  performance  at  the 


40  80  120  160  200 

Current  density  (mA/cm2) 


Fig.  7.  Net  power  output  under  different  operation  conditions  (actuation 
frequency  =  200  Hz). 


fuel  flow  rate  of  5  ml  min-1  is  improved  but  the  net  power 
output  is  still  too  small.  The  other  way  to  achieve  a  high 
flow  rate  with  a  reasonable  net  power  output  is  by  using  two 
micropumps  connected  in  parallel.  The  performance  of  two 
micropumps  with  a  flow  rate  of  2.5  ml  min-1  each  is  also 
shown  in  Figs.  6  and  7,  which  is  much  better  than  that  of 
one  micropump  with  a  flow  rate  of  5  ml  min-1 .  However,  the 
trade-off  is  that  the  material  cost  and  the  system  complexity 
are  increased. 

When  the  driving  frequency  is  equal  to  the  resonance 
frequency,  the  displacement  of  the  diaphragm  will  reach 
the  maximum  value.  After  that,  to  increase  the  driving  fre¬ 
quency  further  will  not  improve  the  pump  performance.  Also 
it  should  be  noted  that  during  the  high  frequency  range  (close 
to  the  resonance  frequency),  our  quasi-static  model  may  not 
be  applicable  and  a  dynamic  model  need  to  be  established. 

Eq.  (31)  together  with  Eq.  (16)  also  describes  the  effects 
of  other  factors  on  the  pump  power  consumption.  All  these 
factors  are  related  to  the  dimensions  and  material  properties 
of  the  pump  diaphragm  structure.  It  is  easily  shown  that  larger 
R ,  ds  i  and  ij  values  are  preferred.  Smaller  dielectric  permittiv- 
ity  ds  of  piezoelectric  material  corresponds  to  smaller  power 
consumption.  However,  the  influence  of  thickness  is  other¬ 
wise  quite  complicated.  Generally  speaking,  for  a  specific 
piezoelectric  disk  thickness,  an  optimal  diaphragm  thickness 
exists.  For  some  detailed  study  on  the  dimensions  and  ma¬ 
terial  properties  of  the  pump  actuator,  please  refer  to  Li  and 
Chen  [12]. 

6.  Conclusions 

In  this  paper,  a  viable  fuel  delivery  system  using  a 
valveless  piezoelectric  micropump  in  a  miniaturized  direct 
methanol  fuel  cell  (DMFC)  device  is  presented.  A  simple  an¬ 
alytical  pump  model  combined  with  some  experimental  data 
from  the  DMFC  device  is  used  to  evaluate  the  performance  of 
the  fuel  cell  system.  Based  on  theoretical  analysis,  the  effect 
of  piezoelectric  materials  properties,  driving  voltage,  driving 
frequency,  nozzle/diffuser  dimension  and  other  factors  on  the 
performance  of  the  fuel  cell  system  is  discussed.  It  has  been 
demonstrated  that  with  an  optimal  system  design,  the  piezo¬ 
electric  micropump  only  consumes  a  very  small  amount  of 
power  generated  by  the  fuel  cell  device.  A  large  driving  fre¬ 
quency  is  suggested  to  decrease  the  power  consumption  of 
micropump.  However,  both  model  improvements  and  exper¬ 
imental  work  are  needed  for  further  study  in  the  future.  It  is 
expected  that  the  design  combining  DMFC  device  with  piezo¬ 
electric  micropump  for  fuel  delivery  leads  to  a  compact  struc¬ 
ture  of  the  power  generation  system  for  portable  electronics. 
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